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Quark coalescence predictions for various multi-strange baryon to anti-baryon ra- 
tios for central Au+Au collisions at RHIC energies are compared to preliminary 
data of the STAR collaboration. The formation of Quark Matter and the sudden 
recombination of its constituent quarks into hadrons is found to be in agreement 
with these preliminary data. It seems that strange hadron and antihadron produc- 
tion in Au-|-Au collisions at RHIC is similar to that in Pb+Pb collisions at CERN 
SPS. 



In the Quark Matter'99 conference the npssibihty of an interesting sce- 
nario for the RHIC reactions was pointed out .El Namely we may meet the situ- 
ation that in the first stage of the colhsion a sort of quark gluon plasma fireball 
is formed. This fireball expands and cools, and just before the hadronization 
this matter may change into a Quark Matter similar to that state which was 
present at the SPS Pb-|-Pb collisions at the beginning of hadronization. The 
basic degrees of freedom of this state are the dressed valence quarks, the gluz 
onic degrees of freedom are suppressed e.g. due to their large effective massfl 
The earlier fragmentation of these gluons resulted in new valence quarks, 
that are the constituents of the final state hadrons. Keeping in mind this 
possibility, theoretical predictions were made for the RHIC Au+Au collisions 
based on a coalescence model, which succesfuUy described - among others - 
the antihyperon to hyperon ratios at SPS. 

Recently prelimiaai-y experimental data were shown for these ratios by the 
STAR collaboration.cla Thus it is timely to compare these new experimental 
results with the theoretical predictions. 

1 The Linear Coalescence Model 

In the following we shall demonstrate the deficiency of the linear coales- 
cence model in a very simple case. The idea of coalescence was created many 
years ago for the understanding of some low energy nuclear physics problems. 
Among others the production of deuterons in low energy heavy ion reaction 
was explained by the assumption that a neutron and a proton will coalesce if 
they are near to each other in space and momentum space. We shall denote 
by k the number of particles of type k within the rapidity interval, in which 
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the particle ratios are constant. (In low energy reactions these numbers are 
the total numbers of particles within the fireball.) The number of produced 
deuterons, d, is assumed to be proportional to the number of protons, p, and 
the number of neutrons, n. The factor of proportionality, Cd, can be calcu- 
lated from the elementary cross sections and the phase-space distribution of 
protons and neutrons. This leads to the linear coalescence model prediction 

d — Cdpn. (1) 

This kind of relationship expresses the essence of the linear coalescence model. 
It is an acceptable approximation as long as a small fraction of protons and 
neutrons will coalesce into deuterons, Cd << 1. The linear coalescence model, 
however, has a serious deficiency if a substantial part of the nucleons are forced 
into deuterons. To demonstrate the problem, let us assume, that p — n = A, 
and all the nucleons have to be confined into deuterons at the end of the 
fireball lifetime. This way we model a situation similar to the hadronization 
of the constituent quark matter. The conservation of baryon number demands 
that the baryon charge should be preserved during the coalescence process: 

p + n = 2d (2) 

which leads to the requirement 

p + n — 2cdpn, (3) 

that implies 

2A = 2cdA^. (4) 

Baryon number conservation and linear coalescence thus requires 

Cd - l/A. (5) 

However, the coalescence coefficient Cd can be calculated from the cross section 
of the p + n —^ d reaction and the phase-space distribution of protons and 
neutrons, independently of the total number of these particles, 2A. Hence 
eq. cannot be fuUfiUed in the general case and the assumption of linear 
coalescence leads to contradiction. 

What is the logical error in the linear coalescence model which leads to 
the above, obviously unacceptable result? Essentially this is the lack of a 
feed-back: those nucleons, who are already bound in a deuteron, cannot be 
used any further for the formation of new deuterons. Out of A protons and 
A neutrons, exactly A deuterons can be formed, if a bound state formation 
happens with unit probability. However, linear coalescence calculation would 
predict cx A^ deuterons, which is way too large, if A >> 1. 
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This sort of error in linear coalescence can be corrected by introducing a 
suitable renormalization factor, 6, which takes care of the additional reduction 
of produced particles due to the lack in the reactivity of the already bound 
constituents. This leads to reconciliation of the constraints from conservation 
laws and the proportionality of produced particles with the available number 
of their constituents. 

With this re-normalization factor, the coalescence equation eq. (|l|) is mod- 
ified, 

d = Cd bpp bnU , (6) 

and eq.(Q) still ensures the conservation of the baryon charge. In the simple 
case of p = n = A, this determines the normalization factor, bp = = b and 
predicts the of deuterons, d, as follows: 

A^pn = A^. (8) 

Note also that the total number of deuterons after the completion of the 
recombination becomes independent from the exact value of the coalescence 
coefficient Cd, because of the constraint that all protons and neutrons have 
to be converted to deuterons in this example. This has to be contrasted to 
the result of linear coalescence calculations, where the formation of a small 
number of deuterons from a gas of protons and neutrons depends essentially 
on the coefficient of coalescence, Cd. 

The method of introduction of normalization coefficients can be gener- 
alized for more complex cases that involve the full recombination of differ- 
ent type of constituents into various kind of bound states subjected to con- 
servation laws. In case of hadronization (recombination of valence quarks 
into hadronsl,_tlie resulting non-linear quark coalescence model is the AL- 
COR modelpBfl the name of which stands for ALgebraic COalescence in 
Rehadronization. 

2 ALCOR: Quark Combinatorics in Rehadronization 

The nonlinear ALCOR coalescence model was created for situations, where 
the subprocesses are not independent, they compete with each other. In 
this model the coalescence equations relating the number of a given type of 
hadron to the product of the numbers of different quarks from which the 
hadron consists reads as: 
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p ^ Cpbqbqbqqqq 
A = Ca 5q 6q 6s g g s 
E^C'Ebqbsbsqss 
Q = Cnbs bs bs s s s 



P = Cpb^bqbqqqq 
A = Cj^bqb^bsqqs 
E — C^b^bsbsqss 
Q = CTrbsbsbsSss 



TT = C^bqbqqq 
K ^Ck bqbsqs 
K = C-j^bqbsqs 

7] = Cr,bsbsSS (10) 

Here the normalization coefficients, bi, are determined uniquely by the 
requirement, that the number of the constituent quarks do not change during 
the hadronization — which is the basic assumption for all quark counting 
methods: 

s = Sn + 2E + A + l{ + r] 
s = 3Ti+2E + A + K + 1^ 
q = 3p + E + 2A + K + n 

q^3p + E + 2A + K + n . (11) 

In eq. (^l|) TT is the number of directly produced pions. (Most of the 
observed pions are created in the decay of resonances.) Substituting eqs. ^ 
^ into eq. (0) one obtains equations for the normalization constants. These 
constants are then given in terms of quark numbers and Cj factors. 



3 ALCOR Predictions Compared to Preliminary RHIC Data 

3. 1 Parameter independent hyperon ratios for SPS and RHIC 

In order to arrive at the detailed predictions of the ALCOR model, one has to 
solve numerically the set of equations eqs. (P-pl]). The ALCOR predictions for 
the particle multiplk;ities (especially for the antihyperon to hyperon ratios), 
were shown in RefsJj'H for the SPS energy. 

Later on an interesting observation was pointed out in RefJj Namely it 
was shown, that within the framework of the linear coalescence model in the 
antiparticle/particle ratios the uncertain Cj coalescence factors drop out, and 
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thus one obtains direct relations between particle ratios. Unfortunately the 
linear coalescence model contradiscts to the different conservaticai. laws. 

In the later development, however, it was shown in Refs. Era that even 
in the charge conserving normalized (non-linear) coalescence model one can 
obtain parameter independent ratios between particle multiplicities. These 
relations have the form: 



h-h (-) 

where 

K+ 

D-T^ . (13) 

Thus if the dynamics of hadronization is the coalescence process, then 
these relations should hold among-^the measured quantities. The preliminary 
results of the STAR experiment Uu for Au+Au collision at ^/s — 130 AGeV, 

p/p = 0.61 ± 0.06 

A/A = 0.73 ±0.03 

E/E = 0.82 ±0.08 

i^+Zif- = 1.12 ±0.06 (14) 



clearly satisfy the relations eqs.(|12|-|13D, as required by quark combinatorics. 

Similar results support the quark coalescence mechanism- jabserved in 
Pb+Pb colhsion at i;bcam = 158 AGeV energy at CERN SPsMo 

p/p = 0.07 ±0.01 
A/A = 0.133 ±0.007 
E/E = 0.249 ±0.019 

0/fi = 0.383 ±0.081 (15) 

which also satisfy the relations among the ratios demanded by quark coales- 
cence with 

K+/K- = l.8±0.2. (16) 

Hence the multi-strange baryon to anti-baryon ratios in central Au+Au 
collisions at RHIC and in central Pb+Pb collisions at CERN SPS arc con- 
sistent with the formation and the sudden rehadronization of Quark Matter, 
consisting of valence quarks and anti-quarks. 
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3.2 Detailed ALCOR predictions for RHIC 



The first ALCOR predictions for RHIC were given for Au+Au collision at 
^/i = 200 GeV in Ref.EI The first data sets were taken at the bombarding 
energy of -^/s = 130 A GeV, and the preliminary values for particle ratios 
were reported by the STAR Collaboration in October 2000 at the XXXth In- 
ternational Symposium on Multiparticle Dynamics as well as at the Quark 
Matter in January 2001.0. These data were similar to the ALCOR predictions 
for particle ratios at -^s — 200 AGeV. Here we re-evaluate the ALCOR pre- 
dictions for the lower than expected values of ^/s. We utilize the capability 
of the ALCOR model to work in a reasonably small rapidity bin (Ay « 1) 
where all the available quarks and antiquarks are able to interact to form the 
measured mesons and baryons. The calculation of hadron multiplicities in 
the ALCOR model is based on the number of quarks and antiquarks in the 
rapidity interval under study. These quark and antiquark numbers consist of 
the WM, dd, ss pairs produced in the collision and the u and d quarks of the 
colliding nuclei. We utilize the measured (preliminary) value of the central 
rapidity density of the negatively charged hadronsH 

dA^h-/dy = 264± 18 . (17) 

If we apply the parameter of the strangeness (/s = Nss/ {N^u + N^-^) =ri^-^^) 
and baryon production (aeff = 0.85) determined at CERN SPS energy,El then 
we are very close to make a full calculation of all hadronic ratios. The only 
missing information is the number of stopped nucleons in this rapidity region. 
(In previous calculations at CERN SPS energy we have not used this infor- 
mation, because total hadron numbers were calculated including all stopped 
nucleons.) Choosing 197 newly produced uu and dd pairs (this means 86 
ss pairs) and assuming that 5 % of the total nucleon numbers stopped into 
the central unit of rapidity (which means 8 proton and 12 neutron), we can 
calculate the measured rapidity densities and hadronic ratios for the Au+Au 
collision at ^/s = 130 AGeV in the central rapidity as displayed in Table 1. 

Investigating particle production in the central rapidity we repeat our 
calculation for Pb-|-Pb collision also at CERN SPS energy. The central 
rapidity densifai. of the negatively charged hadrons was measured by NA49 
Collaboration:tj 

dA^h-/dy = 195 ± 15 . (18) 

In the ALCOR calculation we use 121 newly produced uu and dd pairs (this 
means 53 ss pairs) and assume 17 % of the total nucleon numbers to be 
stopped into the central unit of rapidity (which means 28 proton and 42 
neutron). For this case ALCOR model yields dN]^-/dy — 195 in the central 
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ALCOR model 


Preliminary data 


Ref. 


h- 


260 


264 ± 18 


STAR [4] 




0.63 


0.61 ±0.06 


STAR [4] 


A/A 


0.72 


0.73 ±0.03 


STAR [4] 




0.83 


0.82 ±0.08 


STAR [4] 


K+IK" 


1.13 


1.12 ±0.06 


STAR [4] 


K+/'K+ 


0.142 


0.15 ±0.01 


STAR [13] 


K*+/h- 


0.077 


0.065 


STAR [4] 




0.067 


0.060 


STAR [4] 



Table 1. Hadron production in Au+Au collision at y/s = 130 AGeV from the ALCOR 
model and the preliminary experimental data [4,13] 



rapidity and the following baryon/antibaryon ratios: p/p — 0.07, A/A = 
0.125, S/S = 0.230 and U/n = 0.425. The obtained ratio K+fR- = 1.73 
fuUfils also the expectation, see eqs.( p^|l^ . 

The calculated and the measured antibaryon/baryon ratios at SPS and 
RHIC energies are displayed in Figure 1. 



4 Summary: Quark Matter at CERN SPS and at RHIC 

It seems that the gluonic degrees of freedom are liberated in the initiajLstate 
in Au+Au collisions at RHIC, as evident from data on jet quenching |l3'Ej and 
this initial state seems to be different from that of Pb+Pb collisions at CERN 
SPS.t£l As the fireball expands and cools the gluon dominated initial state in 
the y/s = 130 AGeV Au+Au reactions at RHIC decays to a quark matter 
kind of state, where constituent quark and antiquark degrees of freedom are 
dominant. This quark-antiquark matter undergoes a sudden hadronization 
similarly to that at CERN SPS. 

The success of the coalescence model to describe the experimental data at 
RHIC and SPS strongly supports the validity of the basic assumption behind 
the ALCOR model: before the hadronization, the basic degrees of freedom of 
the expanding fireball are massive (constituent) quarks and antiquarks. These 
elements of the deconfined phase coalesce into final state hadrons in a sudden 
process. 
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Figure 1. Particle ratios at mid-rapidity in Pb+Pb collisions at SPS (open circle) and in 
Au+Au collision at RI 
square) are from RefsJ; 



Au+Au collision at RHWr-jp«en square). Preliminary experimental data (filled circle and 



Acknowledgments 



This research has been supported by the two Bolyai Fellowships of the Hun- 
garian Academy of Sciences (T. Cs. and L. P.), by the Hungarian OTKA 
grants T024094, T025579, T026435, T034269, by the US-Hungarian Joint 
Fund MAKA 652/1998 and by the NWO-OTKA grant N025186. 



ismdOO-zimanyi: submitted to World Scientific on February 1, 2008 8 



References 



1. J. Zimanyi, Nucl.Phys. A 661 (1999) 224. 

2. P. Levai and U. Heinz, Phys. Rev. C 57 (1998) 1879. 

3. W. Christie for the STAR Collaboration, Proceedings of XXXth 
ISMD'2000, Tihany, Hungary, World Scientific, Singapore, 2001, Eds. 
T. Csorgo, S. Hegyi and W. Kittel. 

4. J. Harris talk at Quark Matter'2001, to be published in Nucl. Phys. A 
(2001). See also ^ittp:// www.rhic.bnl.gov/qm200T7| . 

5. T.S. Biro, P. Levai, and J. Zimanyi, Phys. Lett. B 347 (1995) 6. 

6. T.S. Biro, P. Levai, and J. Zimanyi, J. Phys. G 23 (1997) 1941. 

7. T.S. Biro, P. Levai, and J. Zimanyi, Phys. Rev. C 59 (1999) 1574. 

8. A. Bialas, Phys. Lett. B 442 (1998) 449^ 

9. J. Zimanyi, T.S. Biro, and P. Levai, \iep-ph/9904Ml , 

J. Zimanyi, T.S. Biro, T. Csorgo, and P. Levai, Phys. Lett. B 472 (2000) 
243. 

10. M. Kaneta et al., NA44 coU. J. Phys. G 23 (1997) 1865. 

11. R. Caliandro et al., WA97 coU. J. Phys. G 25 (1999) 171. 

12. S.V. Afanesiev et al. (NA49 CoU.) Nucl. Phys. A 610 (1996) 188. 

13. N. Xu talk at Quark Matter'2001, to be published in Nucl. Phys. A 
(2001). See also ^ittp:// www.rhic.bnl.gov/qm200T7| . 

14. G. David for the PHENIX Collaboration, Proceedings of XXXth 
ISMD'2000, Tihany, Hungary, World Scientific, Singapore, 2001, Eds. 
T. Csorgo, S. Hegyi and W. Kittel. 

15. P. Levai et al., Proceedings of XXXth ISMD'2000, Tihany, Hungary, 
World Scientific, Singapore, 2001, Eds. T. Csorgo, S. Hegyi and W. 
Kittel. 

16. T. Csorgo, Proceedings Torino 2000, |hep-ph/0011339| , Nucl. Phys. B. 
Proc. Suppl. in press. 



ismdOO- zimanyi: submitted to World Scientific on February 1, 2008 9 



